arXiv:1501.06798vl [quant-ph] 27 Jan 2015 


Demonstration of quantum permutation algorithm with a single 

photon ququart 

Feiran Wang, Yunlong Wang, Ruifeng Liu, Dongxu Chen, Pei ZliangQ Ilong Gao, and Fuli Li 
MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, 

and Department of Applied Physics, 

Xi’an Jiaotong University, Xi’an 710049, China. 

Abstract 

We report an experiment to demonstrate a quantum permutation determining algorithm with 
linear optical system. By employing photon polarization and spatial modes, we realize the quantum 
ququart states and all the essential permutation transformations. Compared with the classical case, 
this work determines the parity of the permutation in only one step of evaluation and displays the 
remarkable speedup of quantum algorithm. This experiment is accomplished in single photon level 
and exhibits strong universality in quantum computation. 
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Introduction 

As quantum counterpart of classical computer, quantum computer reveals incredible ef¬ 
ficiency to execute arithmetic tasks and threatens the security of classical communication. 
Quantum algorithm is the sole of quantum computation, and it shows the amazing power of 
quantum parallelism and utilizes quantum interference to get the expected result, therefore 
it attracts particular concern to develop new quantum algorithm in recent years. The con¬ 
cept of simulating physics progresses with quantum computers was originated in Richard 
Feynman’s observation that computers built from quantum mechanical components would 
be ideally suited to simulating quantum mechanics [I]. Since then, the first efficient quan¬ 
tum algorithm was proposed by Deutsch in 1985 [2] and generalized by Deutsch and Jozsa in 
1987 [3j. Lately, an increasing number of practical programs are presented, such as factor¬ 
ing large integer [4J, Grover’s searching algorithm for database [5] and Simon’s exponential 
acceleration algorithm for the black box problem [6]. What’s more, Harrow et al. in Bristol 
came up with a quantum scheme to decrease the computational complexity of solving linear 
system of equations from O(n) to login), and this is the first quantum algorithm to work 
out the most fundamental problems in engineering science [7]. Some quantum algorithms 
have been demonstrated in different physical systems, such as ion traps mu , supercon¬ 
ducting devices [f~2] fhf] . optical lattices [E3, IT5] , quantum dots PH EH], and linear optics 
PSH23]. Linear optical system is potential to realize quantum computer, due to its good 
scalability, easy-handling and high stability, and it is a good candidate for implementing 
quantum algorithms. 

A new algorithm based on quantum Fourier transformation has been proposed to identify 
the parity of a permutation recently [2T]- The permutation problem can be described as a 
black box problem, which is to determine the parity of the permutation based on the n input 
and output states of black box. For example, in a set of six different possible permutations 
(1,2,3), (2,3,1), (3,1,2), (3,2,1), (2,1,3) and (1,3,2), the first three are positive cyclic or 
even permutations while the last three are negative cyclic or odd permutations. If we treat 
permutation operation as a function /(x) defined on the set x G {1,2,3}, determination 
of parity requires the evaluation of /(x) for at least two different input values of x with 
classical algorithm. 

In this article, we consider four orthonormal states, i.e., |1) = (1,0,0,0) T , \2) = 
(0,1,0,0) T , |3) = (0,0,1,0) T , and |4) = (0,0,0,1) T , and the bijection is /: {1,2, 3,4} —s- 
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{1,2, 3,4}. For the input state (1,2,3,4), there are eight different possible output states. 
These eight states and the corresponding transformations are divided into two categories as 
mentioned above. The positive cyclic permutations and the corresponding unitary transfor¬ 
mations are: 
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and the negative cyclic permutations and the corresponding unitary transformations are: 
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The black box can be considered as a special device which can realize the corresponding 
operation for a given permutation task. 

To determine the parity, a direct running of the operator on the eigen states of ququart 
will not work. Therefore at least twice runnings are needed to evaluate the box parity as 
in the classical case. For example, if we input a state |2) and get the output state |4), 
this progress corresponds to two possible permutation transformations f 3 and f$. So to 
determine the parity of permutation, we need at least another running. 
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FIG. 1 . Quantum circuit to realize the permutation algorithm. Uft denotes the Fourier 
transformation. The black box contains all the eight permutation operators in the algorithm and 
Up T is the inverse Fourier transformation. 
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Fig- [3 gives the quantum circuit to realize the permutation algorithm. Considering the 
quantum algorithm of this task, we start from a superposition state: 


|^n) = |^2) = (|l)+i|2)-|3)-i|4))/2, 
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here | ipj) = Upr\j), (.?' = 1,2, 3,4) which is generated by quantum Fourier transform Uft'- 
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Then the input state get through the black box, and the output states can be written as 
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where k — 1, 2, • • ■ , 8 label eight permutation transformations. Application of U k on the 
input state |^ in ) gives |r/> 2 ) for even f k and it gives |^ 4 ) for odd f k . By employing the inverse 
Fourier transform U FT and checking the final state, we will know the parity of permutation 
is even (odd) when we acquire the state |2) (|4)). Thus, the quantum algorithm allows us 
to determine the parity of a cyclic permutation with single evaluation of the permutation 
function instead of two. 

Reviewing this algorithm, the Fourier transform Uft and inverse Fourier transform U FT 
are not necessary if we can directly prepare the superposition state \ip 2 ) and discriminate 
the output state |^ 2 ) and Luckily, by employing the photon’s polarization and spa¬ 

tial mode, we can easily realize the ququart |^ 2 ) and distinguish |-0 2 ) and |z/q) in our 
experiment. We should address that |^ 4 ) is also a proper input state for the permutation 
discrimination, and in our experiment, we demonstrate the algorithm for both |^ 2 ) and |^ 4 ). 


Results 

In this context, we utilize photon’s polarization and spatial mode to code the ququart ra, 
and we carry out the whole eight different permutations in a black box which is composed 
by Dove prism (DP) and half wave plates (HWP). Combining the polarization and spatial 
mode as composite quantum state, the four states of ququart can be defined as 

\H,I) |1), | H,r) |2>, ICO 3), \V,r) |4>, (6) 
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where H (V) represents horizontal (vertical) polarization and r ( l ) represents right (left) 
spatial mode. The input state |-0 2 ) can be written in the form 

IV’a} — \H,l) +i\H,r) — \V,l) — i\V,r). (7) 


Then the black box will carry out the permutation transformation on |'02)- For example, 
suppose the operation is / 2 , the output state after black box can be expressed as \H,r) + 
i\V,l) — | V,r) — i\H,l), which is equal to — i\ip 2 )- 
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FIG. 2. Experimental setup of permutation algorithm, (a) The single photon source is 
achieved by deeply attenuating coherent light into single photon level. P denotes the polarizer for 
preparing horizontal polarization state. HWP3 are used for initial polarization state preparation. 
Two beam splitters (BS) and four mirrors (M) are used to set up a Mach-Zehnder interferometer. 
The piezo transmitter (PZT) actuator is used to modulate the phase ip between l and r paths. The 
black box is used to realize eight permutation transformations. Two detectors (D1 and D2) are 
single-photon counting module (SPCM-AQRH-14-FC) used to record the count of photons, (b) 
The black box consists of a DP at 0° and two HWPs at 45°. Different permutation transformation 
can be achieve by different combination of DP and HWPs. 

The sketch of the experimental setup is shown in Fig. [2j)a) . The source is achieved by 
deeply attenuating coherent light into single photon level (coherent parameter a ~ 0.1). 
Polarizer (P) and HWP3 are placed to prepare the initial input polarization state and beam 


5 













splitter (BS1) is to prepare spatial state. The black box consists of a DP and two HWPs 
as shown in Fig. [2^b), the dove prism is located at an angle 0° relative to the horizontal 
plane. Photons with different spatial modes r and l pass through DP will swap their spatial 
modes, and the polarization will be altered by HWP at at 45° succinctly. To realize the 
permutation transformations with black box, what we need is adding or removing the DP 
and HWPs. 

In this scheme, we prepare the initial polarization state of photon in a superposition 
state (| H) — |V)) /\/2, then we transform the spatial mode into (| l) + |r)) /x/2 by BS1. The 
relative phase <p between two spatial modes can be adjusted by a PZT mounted on the first 
reflecting mirror (Ml). Thus the input state for the black box can be written as 

Wto) = (\H) ~ \V)) (\l) + e»\r))/2 

( 8 ) 

= (\H, l)+e vp \H, r) - | V, Z) - e‘*| V,r)) /2. 

Especially, when <p = 7t/ 2, the initial state | ipi n ) turns to be l^) which is expressed in 
Eq. (Jt|) ; when (p = 37t/ 2, the initial state \ipin) turns to be \ 1 Jj 4 ). Then the state is injected 
into the black box to carry out the permutation transformations. The states of ^ 2 ) and |^ 4 ) 
in our definition can be discriminated directly by the phase of spatial mode. So different 
spatial states interfere at the second beam splitter (BS2), and the parity of permutation can 
be determined by checking which single photon detector (D1 or D2) clicks. 


TABLE I. Experimental implementation of eight different permutations 
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Implementation of positive cyclic permutation. For the identity operation f\ , we 


6 





















remove the DP and HWPs, the input states pass through the setup remain unchanged and 
we achieve the /i operation. For the f 2 operation, DP and HWP1 are required in the setup. 
When the photons in different paths pass through the DP, they will exchange their paths (l —y 
r and r —y l), and the state | ip in ) changes into | H, r) + e l ‘ p \H, l) — | V,r) —e llf \V, l). After that, 
HWP1 at an angle 45° is located at the l path, it will alter horizontal polarization to vertical 
polarization and vice versa. So far we get the final state | H,r) + e tip \V, l) — \V, r) — e lv \H, l) 
which is equal to —i^) when ip = 7t/ 2 and achieve the mapping operation / 2 . The fe 
operation can be carried out similarly by removing the DP and putting HWP1 and HWP2. 
Analogous to the action above, the last positive cyclic permutation / 4 can be accomplished 
by inserting DP and HWP2. When photons pass through the device for the positive cyclic 
permutation, a simple calculation shows the output state |Vw) = Uk\i>i n ) {k = 1,2, 3,4) as 
follows: 
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Implementation of negative cyclic permutation. Firstly we place both DP and HWPs 
into the optical route to exchange the spatial mode and the polarization. This will be the 
gate. If we remove the DP and only place HWP1, when the photons pass this setup, it will 
go through an / 6 transformation. For the purpose of carrying out f 7 operation, it is easy 
to obtain that only a DP is needed. The last operation can be achieved by only employing 
HWP2 to change the polarization mode on the right route. Similarly, we arrive the final 
output states: 

,«W„ --(BjjB) ■ - (Ujp.) (tft!) 

(M) (aype) ; = (M2) (M). 

( 10 ) 


From the above discussion, we carry out all the eight essential transformations for the 
parity determining algorithm. All the transformations and corresponding implementation 
approaches are summarized in Table |l| In particular, when ip = —, the output states in 
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Eqs. (J9| and (10) are equal to ^ 2 ) and |z/q), respectively. The Eqs. (J9]) and (10) clearly 
show that the polarization of output states are same while the relative phase between two 
spatial modes are different for two types of parity. This means that we can accomplish 
the identify progress of permutation through the spatial bits without the inverse Fourier 
transform. After the spatial modes interfered on BS2, we can determine the parity is odd 
(even) when detector D1 (D2) clicks. With similar analysis, when the input state is | 7 /q) 

r 

(relative phase ip — —), we can determine the parity is odd (even) when detector D2 (Dl) 
clicks. 

We record the photon counts of the two detectors Dl and D2 with 0.5 V a step of the 
PZT voltage synchronously. Our experimental results are shown in Fig. [3] Figure 3(< 


shows the results of positive cyclic permutation transformations from /i to fi, and Fig. 3(b) 


gives the results of the remaining four odd operations. The black square spots represent 
the counts of Dl, and the red triangular spots corresponding to the counts of D2. As we 
discussed above, when the relative phase ip is equal to (2N + 1/2) 7r, where N is an integer, 
only detector Dl clicks for the odd parity and D2 clicks for the even parity. These special 
points are pointed out by the green dashed line, and from these points we can get the parity 
information of the permutation evidently. Analogously, when the relative phase ip is equal 
to (2N + 3/2)7t, we still can determine the parity by the blue dashed line labelled in Fig. 3. 


Discussion 

We define r] = \Cm — Cm\/\CDi + Cm as the contrast ratio to evaluate the accuracy of 
our experiment, where Cm and Cm denote the photon counts of Dl and D2, respectively. 
Theoretically, the contrast ratio tj is equal to 1. In our experiment, The contrast ratio is up 
to 96 ± 2% for all eight cases in Fig. [3j The error and the small shift of the dashed line for 
different permutations mainly comes from the imperfection of optical devices and moving or 
adding DP and HWPs. 

In conclusion, we briefly introduce the quantum permutation algorithm and put forward 
a scheme to implement this algorithm by employing linear optical element. The main point 
to achieve the parity determining algorithm is the bijection operator and recognize the 
parity through the output states from the two detectors. By using composite quantum 
bit to realize ququart, our experiment is greatly simplified both in state preparation and 
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FIG. 3. Experimental data of the parity determining algorithm, (a) positive cyclic permu¬ 
tation operation, (b) negative cyclic permutation operation. Black square dots show the photon 
counts of D1 , and red triangle points show the photon counts of D2. Fitting lines are also shown. 
Green (blue) dashed vertical lines are used to mark the proper points (phases) of the initial states 
(IV^))) which can be used to perfectly discriminate the two kinds of functions. 

state discrimination, where Fourier transformations Uft and U' FT are saved. Although this 
algorithm only provide a two to one speed-up towards classical case, it shows the power of 
quantum parallelism and quantum computation validly, and expand our thinking for more 
efficient algorithms. We noticed that this scheme has also been realized in spin-| NMR 
quadrupolar system with four energy levels mm- 
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Method 


Figure [2] displays the setup needed in the experimental progress. For the specific operation 
such as / 2 , we only need one HWP at 45° in the left route, however, this HWP gives an 
additional phase because of the intrinsic thickness. With the purpose of eliminating this 
extra effects, we put another HWP at 0° in the right path to compensate phase without 
changing the polarization. We also need to calibrate the voltage of PZT to confirm the 
initial state is proper prepared before the experiment proper starts. 
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